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ABSTRACT 

oredlSina8^  S6ekS  b6tter  understandin8  of  the  general  problem  of 
predicting  the  response  of  fixed  underwater  structures  to  a  shock  wave 

an°el!«Mn8n  r°U8h  ^  Water'  Tw°  ldealized  structures  are  considered- 
ricid  rvIrJ  Y're^ra  u6  r±gid  8phere  and  an  elastically-restrained  ’ 

6  Sh°Ck  18  a8Sumed  t0  be  an  exponentially-decaying 
lit  nhr /“J8*  ?  a?  acoustlc  fluld-  Solutions  to  the  exact  equations 

numer?^tn?  analytlcally  for  the  spherical  structure  and  by  Bellman's 
numericai  inversion  procedure  for  both  structures.  Previously  obtained 

Mindf!°-Rif?ri,the  restralned  cylindrical  structure,  simplified  by  the 

cu^es  arp  f1  apProxJm^lon»  are  found  to  be  in  agreement.  Response 

Thl  pff  ^  8lr!I  f°r  b°th  structural  ^Pes  over  a  range  of  parameters. 
The  effects  of  drag  are  discussed. 
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INTRODUCTION 


This  study  continues  the  analysis  that  was  begun  In  reference  1, 
and  seeks  a  better  understanding  of  the  general  problem  of  predicting 
the  response  of  an  underwater  elastically-restrained  structure  to  a 
shock  wave  propagating  through  the  water.  This  problem  Is  of  Interest 
to  those  who  may  want  to  consider  the  design  or  vulnerability  of  fixed 
ocean-floor  structures.  There  are  sufficient  differences  between  this 
problem  and  that  of  a  neutrally  buoyant  submarine,  for  example,  to 
make  one  suspect  that  existing  analytical  and  experimental  data  may  not 
provide  all  the  answers  needed.  The  fixed  ocean-floor  structure  has 
several  unique  considerations:  (1)  It  is  not  free  to  move  with  the 
water  shock  or  the  flow  behind  it*  this  alters  the  loading  and  pro¬ 
duces  a  greater  flow  of  water  past  the  structure  with  accompanying  drag 
forces.  (2)  It  is  not  necessarily  neutrally  buoyant;  it  may  be  denser 
or  less  dense  than  the  volume  of  water  it  displaces.  The  effect  of 
this  relative  density  is  important  in  the  response  of  fixed  structures. 
(3)  It  depends  upon  the  adequacy  of  the  structural  restraining  members. 
The  design  of  these  may  be  a  critical  problem.  (4)  It  is  subjected  to 
a  different  kind  of  absolute  motion  and  acceleration  than  would  be  a 
free  structure  such  as  a  submarine.  The  structure  motions  must  be  pre¬ 
dicted  so  that  criteria  for  stable  positioning  and  shock  isolation  can 
be  met.  (5)  It  is  restrained  by  and  tied  ta  the  ocean  floor.  It  is 
thus  subject  to  loading  caused  by  motion  of  the  ocean  floor  as  well  as 
that  of  the  water  shock. 

Not  all  of  the  above  considerations  can  be  treated  initially.  The 
special  problem  in  the  beginning  stages  of  such  an  analysis  is  the  dif¬ 
ficulty  in  solving  any  but  the  most  idealized  situation.  In  this  study, 
as  in  reference  1,  a  simplified  structure  is  mathematically  modeled,  and 
the  last  of  the  above  considerations  is  not  included.  For  the  water 
shock  an  acoustic  approximation  has  been  used,  which  implies  certain 
assumptions  about  the  shock  strength  and  velocity  as  well  as  the  nature 
of  the  fluid  surrounding  the  structure.  Although  the  acoustic  approxi¬ 
mation  is  known  to  be  inaccurate,  it  is  not  clear  how  large  an  error  it 
contributes.  In  some  cases  the  error  may  be  minor.  If  the  source  of 
the  underwater  disturbance  is  sufficiently  far  from  the  structure,  the 
shock  reaching  the  structure  will  be  traveling  at  acoustic  speed.  Also, 
any  drag  forces  due  to  steady  flow  appear  to  be  quite  small.  For 
stronger  shocks  the  principal  sources  of  error  probably  arise  from  the 
neglect  of  forces  produced  by  non-steady  flow,  turbulence,  and  separa¬ 
tion  of  the  boundary  layer.  It  was  demonstrated  in  reference  1  that  in 
the  case  of  acceleration  the  peak  values  occur  before  the  shock  has 
traversed  the  structure,  so  that  viscous  drag  and  turbulence  have  not 
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yet  developed.  This  would  imply  that  treating  the  water  as  an  acoustic 
fluid  probably  gives  a  good  estimate  of  peak  acceleration.  On  the  other 
hand,  the  peak  velocity  and  peak  displacement  occur  at  later  times  so  it 
might  be  suspected  that  a  larger  error  is  present  in  the  acoustic  pre¬ 
diction  of  these  values. 

The  reason  for  persisting  in  the  acoustic  approximation  is  that  it 
is  tne  most  tractable  analytically.  However,  even  the  acoustical  equa¬ 
tions  become  very  difficult  to  solve  as  the  number  of  parameters  in¬ 
creases.  In  this  study  at  least  two  additional  parameters  are  of  inter¬ 
est:  The  relative  buoyancy  of  the  structure  and  the  restraint  function. 
A  large  number  of  parameters  makes  it  impracticable  to  run  direct  numer¬ 
ical  solutions  of  the  equations.  It  is  more  desirable  to  try  to  obtain 
final  solutions  as  some  explicit  function  of  the  problem  parameters  so 
that  results  can  be  obtained  for  any  set  of  parametric  values.  This  was 
the  approach  followed  by  Murray2  for  the  case  of  an  unrestrained  cylin¬ 
der;  but  it  did  not  lead  to  a  final  analytical  solution,  and  numerical 
inversion  of  Fourier  transforms  was  required. 

In  reference  1  the  authors  considered  the  rigid  body  response  of 
an  elastically-restrained  cylindrical  deep-ocean  structure  to  an  acous¬ 
tic  shock.  For  this  problem  the  method  employed  by  Murray  becomes  very 
tedious  and  an  effort  was  made  to  find  simple  procedures  of  comparable 
accuracy.  Solutions  were  obtained  using  the  Mindlin-Bleich3  approxima¬ 
tion  which  makes  the  equation  simpler  by  assuming  a  plane  wave  character 
for  the  scattered  and  radiated  waves.  For  the  special  case  of  no 
restraint  these  solutions  were  compared  with  Murray's  results  and  found 
to  be  in  close  agreement. 

It  was  not  possible  to  determine  the  accuracy  of  the  Mindlin- 
Bleich  approximation  for  restrained  cylinders  because  no  solutions  have 
been  found  for  reference.  One  of  the  objectives  of  this  study  has  there¬ 
fore  been  to  develop  such  solutions  to  the  exact  equations  for  the 
restrained  cylinder.  This  is  difficult  because  solution  of  the  exact 
equation  leads  to  transformed  equations  that  cannot  be  inverted  analyti¬ 
cally.  It  is  necessary  to  use  a  numerical  inversion  procedure  as  did 
Murray.  His  procedure  could  be  followed,  but  in  this  more  difficult 
problem  that  procedure  becomes  almost  prohibitively  laborious  and  time- 
consuming.  Instead,  Bellman's  numerical  inversion  procedure4  wac  tried. 
The  appropriateness  of  doing  this  was  demonstrated  for  the  case  of  a 
spherical  structure,  for  which  analytical  inversion  of  the  transformed 
equation  can  be  performed  and  a  comparison  of  results  is  possible. 

This  study  therefore  accomplishes  the  following: 

1.  It  presents  an  exact  solution  for  the  rigid-body  response  of 
an  elastically-restrained  spherical  structure  in  an  acoustic  fluid. 

2.  It  demonstrates  the  appropriateness  of  the  Bellman  procedure 
for  inverting  the  Laplace  transforms  in  the  case  of  the  spherical 
structure  by  comparison  with  1,  above. 
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3.  For  the  case  of  a  restrained  cylindrical  structure  in  an 
acoustic  fluid,  it  presents  a  solution  of  the  exact  equations  using 
Bellman's  procedure. 

4.  It  compares  solutions  using  the  Mindlin-Bleich  approximation 
for  the  cylinder  with  the  results  of  3,  above. 

The  following  sections  discuss  these  results  in  detail  and  result¬ 
ing  figures  show  the  calculated  response  of  both  spherical  and  cylin¬ 
drical  structures.  The  Appendices  contain  the  procedures  used  for 
computing  inverse  transforms. 


SPHERICAL  STRUCTURE 

It  is  possible  to  obtain  an  analytical  solution  of  the  equations 
for  a  restrained  spherical  structure.  Once  that  solution  is  at  hand, 
it  can  be  used  as  a  standard  for  testing  numerical  methods  of  solving 
the  same  equations  in  order  to  determine  the  appropriateness  of  numeri¬ 
cal  solutions.  If  a  satisfactory  numerical  procedure  exists  for  the 
case  of  the  restrained  spherical  structure,  it  is  reasonable  to  try  the 
same  numerical  procedure  in  the  case  of  the  restrained  cylindrical 
structure,  for  which  an  exact  analytical  solution  has  not  been  found. 

Suppose  that  an  elastically-restrained,  spherical  deep-ocean 
structure  is  subjected  to  an  underwater  shock  wave.  The  source  of  the 
shock  wave  is  assumed  to  be  sufficiently  far  from  the  structure  so  that 
the  incident  shock  wave  form  is  a  plane  wave  describable  by  linear 
acoustic  theory.  For  simplicity  the  effects  due  to  reflections  from 
the  ocean  floor  and  the  ocean  surface  are  not  considered. 

Figure  1  describes  the  geometry  and  shows  the  coordinate  system. 
The  *op  diagram  is  for  time,  t  ■  0.  The  shock  wave  is  traveling  from 
left  to  right  and  time  is  measured  from  the  instant  the  shock  wave 
reaches  the  sphere.  Ihe  sphere  is  of  radius  a  and  standard  spherical 
coordinates  are  used;  r  is  the  distance  from  the  origin,  and  0  is  the 
angle  measured  from  the  polar  axis,  which  is  taken  to  be  in  the  direc¬ 
tion  of  propagation  of  the  wave.  The  problem  has  asimuthal  symmetry, 
so  that  a  second  angular  coordinate  is  not  required  to  describe  the 
velocity  field.  The  lower  view  in  Figure  1  shows  the  configuration  at 
a  later  time  t.  The  shock  wave  has  advanced  the  distance  ct  (c  is  the 
acoustic  velocity  of  the  fluid)  and  the  sphere  has  moved  its  center 
through  a  distance  z  in  the  propagation  direction.  For  clarity,  the  z 
displacement  has  been  grossly  exaggerated.  Not  shown  on  the  lower 
diagram  are  the  reflected,  diffracted  and  radiated  waves. 

Three  equations  describe  this  dynamical  situation.  The  first 
expresses  equilibrium  of  forces  acting  on  the  sphere.  These  forces  are 
the  inertia  force,  the  elastic  restraining  force,  and  the  forces  due  to 
the  acoustic  pressures  of  the  incident  and  scattered  and  radiated  waves. 
The  second  equation  is  the  wave  equation  that  describes  the  propagation 
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of  the  acoustic  waves.  The  third  equation  imposes  the  condition  that 
the  water  particle  velocity  at  the  surface  of  the  sphere  must  have  a 
zero  normal  component  relative  to  the  sphere.  Physically  this  states 
that  the  surrounding  water  cannot  penetrate  into  the  sphere. 

The  total  pressure  force  in  the  propagation  direction  is  obtained 
by  integrating  the  incremental  force  in  the  z  direction  due  to  the 
pressures  acting  on  a  strip  of  width  d0  at  6.  This  is  shown  in  Figure 
2.  The  area  of  the  strip  at  angle  6  is  simply  2ira2  sin0  d6.  This 
gives  the  pressure  acting  on  the  strip;  the  component  in  the  propagation 
direction  is  obtained  by  projecting  this  pressure  force  in  the  z  direc¬ 
tion.  The  total  horizontal  pressure  force,  therefore,  is 

Fp(a,0,t)  *  -  ^2ira2  sin0  cos0  d0jp(a,0,t)  (1) 


where  p(a,0,t)  is  the  pressure  distribution  at  the  surface  of  the  sphere 
at  time  t.  The  pressure  itself  is  described  through  the  velocity  poten¬ 
tial  $  for  the  particle  velocity  distribution.  In  terms  of  <b  the  water 
velocity  and  pressure  are  given  by 


Wr(r,9,t)  =  -  —  <|>(r,0,t)  (2) 

p(r,0,t)  =  u  (r ,0 ,t)  (3) 


where  y  is  the  density  of  water.  The  acoustic  field  is  divided  into 
the  incident  acoustic  field  plus  the  scattered  and  radiated  waves.  The 
incident  pulse  is  assumed  to  have  peak  pressure  pD  and  to  decay  with 
time  constant  q  as  given  below. 


-  f2  +  a  -  tl 

P1(r,0,t)  -  pQ  e  L  c  J 

A  „(t .  JsiiJ) 

(4) 

where 

z  =  r  cos0 

In  terms  of  this  incident  wave,  the  equations  for  the  system  are: 


m  2(t)  +  k  z(t)  ■  -  2TTa2|  sin0  cos0 


[Pl(a,e.t)  +  y|l(a*e  »t)Jd6 


V2<ji(r,0,t)  -  77 


1  324(r,0,t; 


i(t)  cos0  -  ^i(r‘0>t)  ii(rfe,t) 

Z(t)  CO80 - 37  -  at 

r  *  a  r  =  a 


The  above  equations  can  be  made  dimensionless  by  defining  the  following 
parameters : 


(£)• 

f  ■V') 


l<j>(r,e,t) 


*-(t)  i 

%>-m2 


s  _  Vs /  _  density  of  structure 

y  density  of  water 


With  these  transformations,  the  three  equations  for  the  system  become: 
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+0)02  ?(t)  =  -  Is  I  sin6  cose  d6  £h(t  -  1  -  cos6)e  "  X(l  "  1  ”  CO80) 

~  0 


+  (i,e,T) 

3t 


(8) 


s2*i „  i  !li  *  i  92*i 

3p^~  +  P  3p  p7  302 


(9) 


dS(t) 

dt 


COS0 


H(t  -  1  -  CO80 ) COS0 


-  X(t  -  1  -  cos0) 
e 


1 


(p»8»t) 


P 


1 


(10) 


Rather  than  deal  directly  with  the  time-dependent  functions,  it  is 
convenient  to  consider  the  Fourier  transforms  of  £(t)  and  The 

general  relations  connecting  a  function  f  (t)  and  its  Fourier  transform 
F  (to )  are : 


f  (t) 


4 


e^01"  F(ui)  du> 


F  (uj  ) 


4 


e  ia)t  f(t)  dt 


Using  this  transformation,  the  above  set  of  equations  become: 
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. .. 


(o)Q2  -  u>2)  £(aj)  “  -  —  I  sine  cose  d0  i= 


-  ±oj t (1  +  cos0) 
X  +  Jw 


+  iu)^1(i,e,oj)  (ii) 


a2  X  i  3^ ,  (p  ,  0  ,to )  1  32ij)  (p,6,0)) 

-±  (p>e>u)  +  - - - - +  _ — - - «  -  u  ♦jCp.e.o.)  (12) 


p2  302 


ib>  C(w)  cos  6 


-  lw(l  +  CO80 )  a  «  x  /  a  \ 

e _ cos0  _  3  4>  i  (P  » 6  »w  ) 

X  +  iu  3p 


P  =  1 


where  c’(oj)  is  the  Fourier  transform  of  £(t)  and  <{>^,0,0))  is  the  Fourier 
transform  of  <J>(p»0»t).  The  middle  equation  is  now  in  the  form  of  the 
Helmholtz  Equation  in  spherical  coordinates  with  asimuthal  symmetry. 

The  Eigen  functions  of  this  equation,  which  represent  outgoing  spherical 
waves,  are 

h£^  (wp)  P ^ (cos0 ) 


where 


h^^(a)p)  *  J,1"*1  order  spherical  Hankel  Function 
of  the  second  kind 

P^(cos0)  ■  order  Legendre  polynomial 


The  general  solution  for  the  velocity  potential,  therefore,  is  given  by 
a  linear  combination  of  these  Eigen  functions,  i.e. 


$j(p 


oo 

s 

a  ■  o 


A£(u)  h^  ;(up)  P^(cos0) 
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The  coupled  sets  of  equations  involving  T(w)  and  A^Cto)  are  obtained  by 
inserting  the  above  expression  into  the  remaining  two  equations  above, 
multiplying  by  P£(cos0)  and  integrating  with  respect  to  6  from  zero  to 
it.  £(to)  appears  only  in  the  l  =  1  equation,  which  corresponds  to  the 
rigid-body  motion  of  the  sphere.  The  higher  angular  modes  would  contri¬ 
bute  to  the  various  deformation  modes  of  the  sphere,  which  are  not  con¬ 
sidered  in  this  report.  The  two  l  -  1  equations  that  result  are: 


(u)0  2  -  to2)  ?(jo)  =  - 


-3±  -  ito 

1  e 


(2) 


2^  xTT^  j(w)  +  ^Ai<w>hi  <“> 


(15) 


ito 


3  e  lwj11(u)) 

e(u,)  s  fir  x'+H)  - 


d[h1(2)(to)] 


dto 


(16) 


The  problem  has  now  been  reduced  to  two  linear  algebraic  equations  in 
two  unknowns,  ^(oi)  and  A^(to).  The  final  expression  for  "c(to)  is: 


?(«) 


-  ito 

e 


T3 

CM 

3 

h1(2)(o>) 

(X  +  ito) 

(oq2  -  to2  +  hj^Oo) 

dto 

d 

h1(2)(to) 

do 

(17) 


This  is  not  quite  as  formidable  an  expression  as  it  appears  because  the 
spherical  Bessel  and  Hankel  functions  are  simple  linear  combinations  of 
senes  and  cosines.  Explicitly  these  functions  are:5 


s  into  coso) 

- n—  -  - 

to*1  (0 


h 


1 


(2) 


(w) 


ito 


The  final  expression  for  ?(to)  is  obtained  by  substituting  these  values 
for  h^(2)(to)  and  j^(to). 
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This  inversion  can  be  done  by  contour  integration  in  the  complex 
u  plane  to  yield  general  solutions  for  certain  cases,  and  numerical 
solutions  for  other  cases.  This  is  discussed  in  Appendix  A. 

Typical  examples  of  the  time-histories  are  presented  in  Figures 
3  through  20.  The  displacement  £,  velocity,  dC/(jT»  and  acceleration. 


,  are  plotted  versus  time. 


T 


for  a  variety  of  the  parameters. 


The  peak  acceleration  does  not  appreciably  change  throughout  the 
range  of  X  and  studied.  As  increases,  the  maximum  acceleration 
in  the  initial  positive  peak  decreases.  However,  the  magnitude  of  the 
first  minimum  increases.  The  net  effect  is  an  almost  constant  peak 
acceleration. 

As  was  to  be  expected,  the  peak  response  increases  as  S  decreases. 

As  was  observed  for  the  unrestrained  case,  the  rate  of  damping  also 
increases  as  S  decreases.  A  denser  structure  experiences  smaller  peak 
responses.  However,  its  motion  will  persist  for  much  longer  times. 

If  the  transformed  equation  for  T(u)  could  not  be  inverted  readily, 
one  could  make  use  of  a  numerical  inversion  procedure.  One  such  pro¬ 
cedure  has  been  developed  by  Bellman  et  al4  for  Laplace  transforms. 

Since  a  Fourier  transform  can  be  restated  to  a  Laplace  transform  by  an 
appropriate  substitution,  the  Bellman  procedure  can  be  applied  to  the 
solution  for  the  sphere.  The  reason  for  doing  this  is  to  obtain  a  check 
on  the  applicability  of  the  Bellman  procedure  to  this  problem  by  compar¬ 
ing  the  results  with  those  obtained  by  contour  integration.  Once  its 
applicability  is  established,  the  Bellman  method  can  be  used  for  the 
solution  for  the  cylinder,  a  solution  that  cannot  readily  be  obtained 
by  contour  integration. 

The  Bellman  procedure  is  described  in  Appendix  B.  This  procedure 
was  used  to  invert  Equation  18,  and  the  results  are  shown  in  Figures 
21  through  30,  which  compare  these  results  with  those  obtained  by  contour 
integration. 
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The  points  shown  in  the  figures  obtained  by  the  Bellman  procedure 
used  a  quadrature  order  N  of  5,  6  and  7,  and  scale  factors  SF  of  1.0 
and  3.33.  The  symbols  used  on  the  plots  are: 


N 

SF 

1.0 

3.33 

5 

X 

+ 

6 

s 

★ 

7 

0 

Y 

The  agreement  for  small  w0  is  very  good.  With  increasing  natural 
frequency,  the  approximation  breaks  down  after  the  first  maximum  is 
reached.  The  main  deficiency  from  a  design  point  of  view  is  that  the 
first  minimum  appears  underestimated  at  these  higher  frequencies.  It 
is  expected,  however,  that  actual  structures  built  on  the  ocean  floor 
will  not  be  that  rigidly  constrained;  hence  aiQ  will  be  relatively  small. 
Bellman's  approximate  inversion  method  should  then  yield  reliable  design 
information. 


CYLINDRICAL  STRUCTURE 

The  next  problem  of  interest  is  the  motion  of  an  elastically  re¬ 
strained  cylindrical  structure  subjected  to  an  underwater  shock  wave. 

As  in  the  previous  section,  the  effects  due  to  the  sea  surface  and 
bottom  are  neglected  and  the  shock  wave  is  considered  to  be  a  plane 
wave  that  is  weak  enough  to  be  describable  by  linear  acoustic  theory. 

The  formulation  for  the  cylindrical  structure  proceeds  in  a  man¬ 
ner  analogous  to  that  for  the  spherical  structure.  The  main  differ¬ 
ences  between  the  two  problems  are  that  the  surface  area  for  the 
cylinder  does  not  contain  the  factor  of  sin9  and  that  the  independent 
variables  for  the  acoustic  fluid  are  now  cylindrical  coordinates  rather 
than  spherical  coordinates.  As  a  consequence  of  these  differences,  the 
Fourier  transform  of  the  cylinder  motion  contains  ordinary  Hankel  func¬ 
tions  rather  than  spherical  Hankel  functions  and  cannot  be  inverted 
analytically. 

Since  the  main  steps  in  the  mathematical  analysis  are  almost  iden¬ 
tical  with  that  of  the  sphere,  only  a  brief  treatment  of  the  development 
will  be  given.  The  upper  drawing  of  Figure  31  shows  the  incident  shock 
wave  just  arriving  at  the  cylinder.  The  motion  of  the  cylinder  as  well 
as  the  fluid  is  independent  of  the  x  coordinate,  which  is  perpendicular 
to  the  plane  of  the  paper.  The  z  and  y  as  well  as  r  and  6  coordinates 
are  shown  and  are  standard  Cartesian  and  Polar  coordinates,  respectively. 
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The  lower  drawing  of  Figure  31  shows  the  situation  at  a  time  t.  The 
shock  wave  has  moved  a  distance  ct  and  the  cylinder  is  displaced  at  a 
distance  z.  As  with  the  sphere,  the  translational  motion  of  the  cylin¬ 
der  has  been  exaggerated  for  ease  of  presentation.  The  equation  of 
motion  for  the  cylinder  is 


where  k  is  the  spring  constant  per  unit  length  of  cylinder  and  m  is  the 
mass  per  unit  length  of  the  cylinder.  The  terms  on  the  righthand  side 
of  the  equation  are  the  pressure  forces  due  to  the  fluid  (the  notation 
is  exactly  the  same  as  with  the  sphere) .  The  wave  equation  for  the 
motion  of  the  acoustic  wave  and  the  boundary  condition  equation  are 
given  below: 


+  i  H  +  1  34  .  1  32  t(r.  e,  t) 

dr7  r  dr  r7  3G7  c7 dtz 


(20) 


«(t)  cose«-f£(r’  e’ 


+  H(ct  -  a  -  a  cos9)  cos0  e 


(21) 


The  transformation  to  dimensionless  variables  proceeds  exactly  as  it  did 
with  the  sphere.  These  transformations  and  the  resulting  equations  follow: 
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c(t) 

^(p.e.T) 


x  - 


s  - 


X(t  -  1  -  cos 6) 


+  TT  (1’e’T 


(22) 


32*j(p.e,T)  1  1  324>1 

— aP -  +  F5r  +  7rae2~ 


92^ 

3^" 


(23) 


d  S.(T) 

dx 


COS0 


H(x  -  1  -  cos6)cos0  e 


-  X(t 


1  -  COS0) 


3^(P  »®»x) 


3p 


P 


1 


(24) 


Uv.ing  the  same  procedure  as  with  the  sphere,  these  equations  are  expressed 
in  terms  of  Fourier  transforms.  The  solution  of  these  transformed  equa¬ 
tions  is  straight-forward,  but  does  not  lead  to  any  expression  involving 
spherical  Hankel  functions  and  Legendre  Polynomials  as  was  the  case  for 
the  sphere.  Instead,  the  solution  for  the  acoustic  field  is  found  in 
terms  of  ordinary  Hankel  functions  and  cos  n0.  The  final  result  for  the 
Fourier  transform  of  the  cylinder  displacement  is  given  below: 


,  4  1 

c(“> 


-  ico 


(icu  +  X) 


co  d 


dco 


CO  H. 


(2) 


-1 


(to) 


-  0)' 


H1(2)(ai)| 


+  CO 


dco 


(25) 
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The  time  dependence  for  the  displacement  is  the  inverse,  which  is  given 
by  the  usual  Fourier  integral.  Due  to  the  fact  that  ordinary  Hankel 
functions  are  much  more  complicated  than  spherical  Hankel  functions ,  it 
is  not  possible  to  perform  the  integration  analytically.  Numerical  in¬ 
version  can  be  accomplished  by  using  the  Bellman  procedure,  which  was 
demonstrated  in  the  discussion  of  the  sphere  and  is  further  explained 
in  Appendix  B. 

Several  examples  for  the  response  of  cylindrical  structures  are 
presented  in  Figures  32  through  43.  The  symbols  have  the  meanings 
previously  defined.  The  displacement,  velocity,  acceleration,  and  time 
are  the  dimensionless  variables  pertinent  to  the  cylinder. 

It  can  be  seen  from  these  figures  that  the  response  of  the  cylinder 
is  similar  to  that  for  the  sphere.  It  is  therefore  inferred  that  the 
reliability  of  the  Bellman  method  extends  to  the  cylindrical  problem  in 
roughly  the  same  parametric  range. 

Fairly  reliable  curves  can  be  drawn  fhrough  these  points ,  and  ad¬ 
ditional  points  can  be  readily  calculated  using  larger  values  of  N  in 
the  Bellman  procedure.  The  values  chosen,  however,  appear  adequate  for 
the  determination  of  the  peak  response  in  the  parametric  range  of 
interest. 

In  reference  1  the  solution  to  the  cylinder  problem  was  obtained 
using  the  Mindlin  and  Bleich  approximation.  Good  agreement  was  obtained 
between  this  approximate  solution  and  Hurray's  exact  solution  for  an 
unrestrained  cylinder.  In  Figures  44  through  46  a  comparison  of  these 
previous  results  with  the  present  calculations  shows  that  the  Mindlin 
and  Bleich  approximation  must  be  reasonably  accurate  for  the  restrained 
cylinder  too.  The  Mindlin-Bleich  results  follow  those  from  the  exact 
calculation  until  the  velocity  becomes  negative.  The  reason  for  this 
is  the  different  treatment  of  the  radiative  force.  In  the  Mindlin- 
Bleich  approximation  this  force  is  always  180  degrees  out  of  phase  with 
the  velocity;  in  the  exact  formulation  the  radiative  force  is  an  integ¬ 
ral  operator  and  bears  no  such  constant  phase  relation  with  the  velocity. 


DRAG  EFFECTS 

Thus  far  the  real  nature  of  water,  namely  its  wetness,  has  been 
neglected.  Due  to  viscosity  the  actual  flow  will  differ  from  that  ob¬ 
tained  in  the  preceding  analysis.  One  immediate  result  is  that  drag 
forces  will  be  present. 

The  drag  forces  due  to  idealized  steady  flow,  however,  would  be 
quite  small  compared  to  the  main  pressure  forces.  For  example,  consider 
the  usual  expression  for  the  drag  force  on  a  sphere.  This  is  given  by: 

FD  -  I*  Cd  1Ta2u2  (26) 
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where  CD  is  the  drag  coefficient  and  U  is  the  relative  speed  of  the 
sphere  through  the  fluid.  Assuming  that  U  ■  UQ  ■  Po/(yc)  *  the 

fluid  particle  velocity  at  the  front  of  the  shock  wave/  the  drag  force 
becomes 


F 


D 


2 


ira* 


(27) 


From  this  one  readily  calculates  the  dimensionless  acceleration, 


(28) 


where 


[41 

F*J» 


is  the  acceleration  due 


be  compared  with  the  main 
example.  For  S  ■  1,  Cn  ■ 


pressure 
2  and  p0 


solely  to  the  drag  term.  This  can 

contribution  for  a  specific  numerical 
-  1,000  psi,  one  gets 


2.25  x  10 


The  calculated  peak  acceleration  from  Figure  32  is  0.7;  therefore,  the 
drag  force  is  about  0.3%  of  the  direct  force.  For  larger  incident 
pressures  this  percentage  increases;  however,  the  linear-acoustic 
theory  itself  collapses  before  this  drag  contribution  is  appreciable. 

In  addition  to  this  small  perturbation,  the  change  in  the  flow 
resulting  from  the  separation  of  the  boundary  layer  may  appreciably 
alter  the  scattering  and  radiative  properties  of  the  sphere.  If  this 
is  the  case,  then  inviscid  theory  could  conceivably  be  in  error  by  much 
more  than  the  0.3%  estimated  above  for  the  contribution  of  direct  drag. 

At  the  present  time  there  seems  to  be  little  information  on  this  subject. 


FINDINGS  AND  CONCLUSIONS 

The  response  of  spherical  and  cylindrical  structures  due  to  a 
weak,  plane  acoustic  shock  wave  has  been  presented.  Exact  expressions 
for  the  Fourier  transform  were  obtained.  For  the  sphere  the  inverse 
was  obtained  exactly;  however,  in  order  to  invert  the  transform  for  the 
cylinder,  an  approximation  was  necessary.  This  approximate  method 
applied  to  a  sphere  was  compared  with  the  exact  solution  for  the  sphere 
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and  found  to  be  acceptable  for  the  range  of  parameters  that  will  be 
encountered  for  a  deep  ocean  structure.  It  can  be  concluded  from  the 
similarity  of  results  for  the  sphere  and  cylinder  that  the  approximate 
method  is  appropriate  for  the  cylinder  in  the  same  parametric  range  as 
for  the  sphere. 

The  direct  drag  effects  due  to  the  viscosity  of  the  water  were 
estimated  to  be  negligible.  Indirect  effects  due  to  the  perturbed 
flow  could  not  be  estimated,  and  a  theoretical  analysis  of  this  effect 
would  be  very  difficult.  A  series  of  small-scale  experiments  should  be 
performed  in  order  to  investigate  the  magnitude  of  this  effect. 


Figure  3.  Displacement  Versus  Time  for  Unrestrained  Sphere;  Step 
Pulse  Incident. 
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Figure  6.  Displacement  Versus  Time  for  Unrestrained  Sphere;  Short 
(X  =  1)  Pulse  Incident. 


Figure  7.  Velocity  Versus  Time  for  Unrestrained  Sphere;  Short  (x  =  1) 
Pulse  Incident. 
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Figure  9.  Displacement  Versus  Time  for  Restrained  (cu  =  2)  Sphere; 
Step  Pulse  Incident.  0 
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Figure  11. 


Acceleration  Versus  Time  for  Restrained  (uj  =  2)  Sphere; 
Step  Pulse  Incident. 
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Figure  15.  Displacement  Versus  Time  for  Restrained  (u) 
Step  Pulse  Incident.  c 
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*  4)  Sphere 
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Figure  16.  Velocity  Versus  Time  for  Restrained  (u)  = 
Step  Pulse  Incident.  0 


8.0 


Sphere; 
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Figure  18.  Displacement  Versus  Time  for  Restrained  (u>  =  4)  Sphere; 
Short  Pulse  (X  =  1)  Incident.  0 
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Figure  20.  Acceleration  Versus  Time  for  Restrained  (u>  =  4)  Sphere; 
Short  Pulse  (x  =  1)  Incident.  0 
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Figure  21.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Unrestrained,  Neutrally  Buoyant  Sphere;  Step  Incident 
Pulse . 
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Comparison  of  Bellman  s  Method  With  Exact  Solution  for 
Restrained  (u)  =0.3),  Neutrally  Buoyant  Sphere;  Step 

Pulse  Incident. 


Figure  23.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Restrained  (<ju  =  1.0),  Neutrally  Buoyant  Sphere;  Step 
Pulse  Incident. 
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Figure  24.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Restrained  (u)  *  1.5),  Neutrally  Buoyant  Sphere;  Step 
Pulse  Incident. 
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Figure  25.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Restrained  (u)  -  2.0),  Neutrally  Buoyant  Sphere;  Step 
Pulse  Incident. 
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Figure  26.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Unrestrained,  Neutrally  Buoyant  Sphere;  Short  Pulse 
(X  *  1)  Incident. 
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Figure  27.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Restrained  (uu  =0.3)  Neutrally  Buoyant  Sphere;  Short 
Pulse  (x  =  1)° Incident. 
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Figure  28.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Restrained  (to  =  1.0)  Neutrally  Buoyant  Sphere;  Short 
Pulse  (x  =  1)° Incident. 


Figure  30.  Comparison  of  Bellman's  Method  With  Exact  Solution  for 
Restrained  ( uo  =  2.0)  Neutrally  Buoyant  Sphere;  Short 
Pulse  (x  =  l)°Incident. 
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Figure  35.  Response  Versus  Time  for  Neutrally  Buoyant,  Restrained 

(uu  =  0.3)  Cylinder;  Step  Pulse  Incident, 
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37.  Response  Versus  Time  for  Neutrally  Buoyant,  Restrained 

(w  m  0.5)  Cylinder;  Step  Pulse  Incident, 
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Figure  39.  Response  Versus  Time  for  Neutrally  Buoyant,  Unrestrained 
Cylinder;  Short  Pulse  (X  =  1)  Incident. 
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Present  solution 


4  Mindlin-Bleich 


Figure  44.  Comparison  of  Mindlin-Bleich  Approximation  With  Present 
Calculation  for  a  Buoyant  (S  =  0.5),  Restrained  (uu  =  0 
Cylinder;  Short  Pulse  (X  =  1)  Incident.  0 


Present  solution  A  Mindlin-Bleich 


Figure  45.  Comparison  of  Mindlin-Bleich  Approximation  With  Present 

Calculation  for  a  Neutrally  Buoyant,  Restrained 

(u;  =  0.33)  Cylinder;  Short  Pulse  (X  =  1)  Incident, 

o 
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Figure  46.  Comparison  of  Mindlin-Bleich  Approximation  With  Present 
Calculation  for  a  Heavy  (S  =  2.0),  Restrained  (u)  =  0.23) 
Cylinder;  Short  Pulse  (X  =  1)  Incident.  0 
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Appendix  A 


INVERSION  OF  FOURIER-TRANSFORMED 
EQUATIONS  FOR  SPHERICAL  STRUCTURE 


The  Integral  for  c(t)  that  appears  in  Equation  18  on  page  9  can 
be  done  by  contour  integration  in  the  complex  to  plane.  In  general  there 
will  be  contributions  from  five  poles,  one  at  to  =  iX  and  the  others  at 
the  roots  of  the  quartic  equation 

ito1*  +  (2  +  to^-i(2+to  2  +  -g-)  u) 2  —  2  u>  ^to  +  2  iwQ2  =  0  (A-l) 


The  roots  of  this  equation  can  be  shown  to  lie  in  the  upper  half  plane, 
as  they  must  to  insure  causality.  It  is  not  possible  to  solve  the 
quartic  equation  in  closed  form  for  arbitrary  S  and  to0  .  However,  for 
a)0  =  0  the  quartic  equation  breaks  down  into  a  double  root  at  to  =  0  and 
two  roots  that  are  the  solutions  of  the  remaining  quadratic  equation. 
For  Uq  not  equal  to  0,  recourse  must  be  made  to  numerical  methods  in 
order  to  obtain  a  particular  solution.  In  terms  of  these  roots,  x^, 
C(t)  is 


Y  _ 

1  -  1  J  (*i  -  Xj) 
j  T  1 
j  - 1 


(A-2) 


where  Xr  =  iq  and  x^,  i  =  1,  ...,  4  are  the  roots  of  the  quartic.  It 
should  be  noted,  that  in  case  of  multiple  roots  this  solution  is  not 
valid. 

It  is  instructive  to  look  at  the  inversion  integral  for  the  case 
co0  =  0  in  some  detail.  For  this  case,  the  above  inverse  of  Equation  18 
becomes 


C(t)  =  Vjr 


-JL  [  _!* 

^1TS  00  W(W 


-  1X)["  -  (>  +  is)1  +  is)  l1  -  ¥)] 


[u  -  (i  +  ^)i  +\/(i  +  is)(  1  -  is) 


(A- 3) 
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Figure  47  shows  the  location  of  the  poles  in  the  complex  go  plane,  as 
well  as  the  two  contours,  Cl  and  C2.  For  negative  values  of  t  the 
exponential  along  the  contour  C2  goes  to  zero.  There  are  no  poles 
within  this  contour;  hence  the  result  is  zero.  As  stated  above,  this 
is  simply  clue  to  causality.  Namely,  before  the  shock  wave  hits  the 
sphere,  it  shouldn't  move.  For  t  greater  than  zero  the  contour  can  be 
closed  in  the  upper  half  of  the  plane.  This  is  Cl  in  Figure  47,  and 
this  contour  encloses  all  four  poles.  The  integral  along  the  real  axis 
is  displaced  slightly  below  it,  so  that  the  pole  at  u  =  0  is  inside  the 
contour.  The  value  of  the  integral  is  2ni  times  the  sum  of  the  residue. 

The  general  feature  of  the  solution  can  be  readily  determined  with¬ 
out  explicitly  presenting  the  cumbersome  algebra.  If  X  =  0,  there  is  a 
first  order  pole  at  id  =  0.  The  residue  at  to  =  0  times  2ni  then  repre¬ 
sents  the  asymptotic  displacement  of  the  sphere,  i.e. 


C(“) 


3/x 

2S  +  1 


If  X  -  0,  then  the  pole  at  to  =  0  is  of  second  order  so  that 


4 


(j) 


T  -*• 


3t 

2S  +  1 


hence 


d  C(t)  _ „  _J _ 

dt  t  •*  “  2S  +  1 

The  poles  at 


yield  exponentially  damped  contributions;  the  damping  increases  in¬ 
versely  with  S.  It  is  interesting  to  observe  that  even  without  a 
restraint  force  these  contain  an  oscillatory  component  for  S  greater 
than  one-half.  It  is  not  obvious  what  the  physical  cause  is  for  this 
oscillation.  However,  it  is  rather  strongly  damped  out  in  time  and 
would  probably  not  be  measurable.  (Murray2  observes  the  same  type  of 
oscillations  for  the  motion  of  an  unrestrained  cylinder.)  For  S  smal¬ 
ler  than  one-half,  the  motion  no  longer  is  oscillatory  but  is  exponen¬ 
tially  damped.  It  is  interesting  to  observe  that  for  S  equal  to  one 
and  X  4  0,  £(°°)  is  equal  to  —  which  is  simply  the  distance  a  particle 

X 

of  water  in  front  of  the  shock  wave  would  move  during  the  decay  time  of 
the  pulse.  For  X  =  0  on  the  other  hand. 
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+  3/  (2S  +  1) 


is. 

dx 


which  for  S  greater  than  one  is  less  than  the  particle  velocity  in  the 
shock  wave.  For  S  =  1  the  sphere  moves  with  the  shock  wave,  and  for 
S  less  than  1,  it  moves  faster  than  the  particle  velocity  in  the  shock 
wave.  This  is  not  an  unreasonable  result,  since  from  potential  theory 
it  is  well  known  that  an  accelerated  sphere  will  lead  or  lag  the  water 
velocity  depending  on  whether  the  density  is  less  than  or  greater  than 
the  water  density.  Since  with  the  potential  theory  there  is  no  force 
on  the  submerged  object  due  to  uniform  flow,  this  motion  will  persist. 
On  the  other  hand  the  inclusion  of  drag  forces  would  maintain  the 
asymptotic  velocity  of  the  sphere  equal  to  that  of  the  water  particle 
velocities . 
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Appendix  B 


NUMERICAL  INVERSION  OF  TRANSFORMED  EQUATIONS 


Although  an  exact  solution  for  the  motion  of  the  sphere  can  be 
obtained,  this  is  not  possible  for  a  cylindrical  structure  because  the 
solution  to  the  transformed  equation  cannot  be  readily  inverted.  A 
numerical  inversion  of  the  Laplace  transform  for  the  cylinder  response 
can  be  obtained,  however,  using  a  method  due  to  Bellman,  Kalaba  and 
Lockett. **  It  is  expected  that  the  motion  of  the  sphere  and  the  cylin¬ 
der  will  be  qualitatively  similar,  so  that  the  limitations  of  Bellman's 
method  in  this  application  can  be  determined  by  applying  it  to  the 
motion  of  the  sphere. 

The  solution  obtained  for  the  sphere  was  based  on  the  Fourier 

transform,  whereas  Bellman's  method  requires  the  Laplace  transform. 

Due  to  the  fact  that  the  system  is  at  rest  for  negative  times,  the 

Laplace  transform  FT  (o)  is  related  to  the  Fourier  transform  F_(u)  by 

L  r 


FL(a)  =  /2V  Fp(-  io) 


In  order  to  demonstrate  this,  let  F  (co )  be  the  Fourier  transform  of 
f(t);  then 


V“> 


"—a 


-  iut 


f(t)dt 


(B-l) 


Then  since  f(t)  =  o  when  t  is  less  than  zero,  this  becomes 


FF<»> 


if 


-  iut 


f(t)dt 


(B-2) 


The  above  expression  for  u  =  -  io  is 


Ff<-  io)  - 


00 

'*1 


e  0t  f(t)dt 


(B-3) 
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However,  the  Laplace  transform  of  f(t)  is 


Fl<0) 


(B-4) 


hence 


Fl(o)  -  ^  Fp(-  io) 


(B-5) 


as  was  to  be  shown.  With  this  relation  the  Laplace  transform  for  the 
spherical  problem  can  be  readily  obtained  from  the  Fourier  transform 
solution. 

To  describe  Bellman's  method  for  inverting  the  Laplace  transform 
let  F(o)  be  the  Laplace  transform  of  f (t) ,  hence 


F(o) 


f(t)dt 


(B-6) 


For  the  numerical  integration,  limits  of  zero  to  one  are  easier  to 
hankie.  This  is  accomplished  by  the  following  change  of  variable. 


y  “ 


(B-7) 


which  when  inserted  into  (B-6)  gives 


F(o) 


f(-  tny)dy 


(B-8) 


This  integral  can  be  approximated  by  an  N-tem  quadrature  formula 
(Bellman  uses  Gauss-Legend re,  although  others  nay  work  as  well). 
Equation  (B-8)  therefore  becomes 
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F(o) 


N 


a 


w 


1  f(-  &nyn) 


(B-9) 


where  the  yn  and  wn  depend  on  the  particular  quadrature  formula  used. 

This  represents  one  equation  for  the  N  unknowns,  (f(-  £nyn) ,  n  =  1,2,...,N). 
N  equations  are  obtained  by  substituting  N  different  values  for  a.  Taking 
this  set  to  be  a  -  1,2,...,N,  equation  (B-9)  yields  the  following  set. 


F(j) 


N 

V  „n  yni  -  1  £(-  *ny„);  J 
n  =  1 


1,2 


,N 


(B-10) 


In  matrix  form  equation  (B-10)  is 


{F}  =  [A ] { f } 


(B-ll) 


where 


{  F }  = 


f(-  tnyir 


{f  > 


f(-  ^nyN); 


The  usefulness  of  the  method  is  due  to  the  fact  that  [A]  ^  is  problem- 
independent.  The  solution  of  equation  (B-ll)  is 


{f}  =  [A]'1  {F} 


(B-12) 
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Numerical  values  for  f(-  £nyn)  are  readily  obtained  once  [A]  is  known. 
For  Gauss-Legendre  quadrature,  Bellman  has  tabulated  [A]“l  and  -  £nyn 
for  values  of  N  from  3  to  15.  It  should  be  noted  that  [A]  is  ill- 
conditioned,  hence  double  precision  (16  digits)  is  needed  to  carry  out 
the  Indicated  operations  in  equation  (B— 12) . 

Since  t  *  -  Jlny,  N  values  of  f(t)  can  be  obtained  by  solvinp  equa¬ 
tion  (B— 12)  for  N  values  of  f(-  £ny).  Unfortunately,  the  tn  values 
are  rather  small  and  tend  to  cluster  near  t  =  o  for  the  range  of  N 
given  by  Bellman.  This  problem  can  be  avoided,  however,  and  f(t) 
determined  for  larger  values  of  t  by  using  the  fact  that  the  Laplace 
transform  of  f(oit)  is  given  by  —  F(o/a)  where  F(a)  is  the  Laplace 
transform  of  f(t).  Therefore,  equation  (B-12)  becomes 

N 

y]  A^'1  F(j/a)  (B-13) 

Other  methods  to  extend  the  time  domain  are  discussed  in  reference  4 
and  will  not  be  discussed  here. 


{(oti)  .  i 
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LIST  OF  SYMBOLS 


Radius  of  sphere  or  cylinder 
Acoustic  velocity  in  water 
Hankel  function,  order,  2nc*  kind 
Bessel  function 
Summation  index 

"Spring"  constant  of  restraint 
Summation  index 

Unit  mass  of  sphere  or  cylinder 
Summation  index 
Unit  pressure 

Peak  pressure  of  incident  pulse 
Time  decay-constant 
Radial  distance 

Ratio  of  structure  density  to  water  density 
Time 

Water  particle  velocity  and  acceleration 
Water  particle  velocity  at  front  of  shock  wave 
Weighting  factor 

Displacement,  velocity,  acceleration  of  structure 
in  direction  of  shock  propagation 

Frequency  dependent  amplitude  of  £  mode 

Drag  coefficient 

Mass  coefficient 

Drag  force 

Inertia  force 

Heavyside  step  function 

Quadrature  order,  Bellman  method 

Legendre  polynomial,  order 

Scaling  factor 

Spherical  coordinate  measured  in  vertical  plane 


c 

T 

0 

U 

P 

a 

T 

J 

X 

co 

CO 

o 


Dimensionless  translational  distance 
Fourier  transform  of  C 
Spherical  or  cylindrical  coordinate 
Mass  density  of  water 
Dimensionless  radial  distance 
Laplace  variable 
Dimensionless  time 
Potential  function 
Fourier  transform  of  $ 

Dimensionless  time  constant 
Angular  frequency  (rad/sec) 

Fundamental  frequency  of  structure  in  air 
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